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Abstract 

We consider a simple model system supporting stable solitons in 
two dimensions. The system is the parametrically driven damped non- 
linear Schrodinger equation, and the soliton stabilises for sufficiently 
strong damping. The purpose of this note is to elucidate the sta- 
bilisation mechanism; we do this by reducing the partial differential 
equation to a finite-dimensional dynamical system. Our conclusion is 
that the negative feedback loop occurs via the enslaving of the soliton's 
phase, locked to the driver, to its amplitude and width. 

1. When a liquid layer is subjected to vertical vibration, a one- or two- 
dimensional periodic pattern forms on its surface. This phenomenon has been 
known since the celebrated Faraday resonance experiment pQ; more recently, 
it was found that the vertical vibration is also capable of sustaining localised 
2D states. These spatially localised, temporally oscillating structures — com- 
monly referred to as oscillons — were observed on the surface of granular 
materials [2] , Newtonian [31 E] and non- Newtonian [H] fluids. Subsequently, 
stable oscillons were reproduced in numerical simulations within a variety of 
models, including the order-parameter equations [3 E] , discrete-time maps 
with continuous spatial coupling [7], semicontinuum [5] and hydrodynamic 
P theories. Although these simulations accounted for the formation of os- 
cillons in several particular physical settings, they did not uncover the core 
of the mechanism which makes them immune from the nonlinear blow-up 
and dispersive broadening. The fact that stable oscillons occur in diverse 
physical media and in mathematical models of various nature, suggests that 
this mechanism is simple and general. It should operate whenever one has 
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a balance of dispersion and nonlinearity on one hand, and of damping and 
phase-sensitive amplification on the other. 

In order to crystallise the main ingredients of this mechanism, a sim- 
ple model of nonlinear distributed system exhibiting parametric resonance 
was proposed recently ^Hj. The model comprises a two-dimensional lattice 
of diffusively coupled, vertically vibrated, damped pendula. In the present 
note we consider the associated amplitude equation whose stationary soliton 
solutions furnish the slowly varying amplitudes of the lattice oscillons. Un- 
derstanding how these 2D stationary solitons manage to resist the nonlinear 
blow-up or dispersive decay in the amplitude equation will provide insights 
into the stabilisation of oscillons in vibrated media. 

2. The amplitude equation we consider, 

itfj t + V 2 ^ + 2\^\ 2 tfj-tp = h^j* -ijip, (1) 

is the parametrically driven, damped nonlinear Schrodinger (NLS) equa- 
tion. Here V 2 = d 2 /dx 2 + d 2 /dy 2 . Apart from the pendulum lattice, eq.(JTJ) 
serves as an amplitude equation for a wide range of nearly-conservative two- 
dimensional oscillatory systems under parametric forcing. Physically, it was 
used as a phenomenological model of nonlinear Faraday resonance in fluids 
[TTJ[T21I1]. Independently, it appeared in the context of optical parametric 
oscillators [T3| . 

Two stationary radially-symmetric soliton solutions are given by 

^± = A±e- w +n(A ± r), (2) 

where r 2 = x 2 + y 2 ; 

A\ = 1 ± \jh 2 — 7 2 , 9 + = - arcsin Cj^j , 9_ = n — - arcsin Cj^J , 
and lZ{r) is the bell-shaped (monotonically decreasing) solution of equation 

TZ rr + -K r - K + 2K 3 = 0, (3) 

r 

with the boundary conditions TZ r (0) = TZ(oo) = 0. This solution is well 
documented in literature [T3] . 

The soliton ip + exists for all h > 7 while the soliton ^_ exists only in the 
wedge 7 < h < y/1 + 7 2 . It is pertinent to add here that when h < 7, all 
initial conditions are damped to zero. This follows from the rate equation 

d t \ij\ 2 = 2V(H 2 V X ) + 2H 2 (/isin2 X - 7), (4) 
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Figure 1: Stability diagram for the two-dimensional soliton (ip+). No stationary 
localised solutions exist below the straight line h = 7, while above the 'parabola' 
h = a/1 + 7 2 all localised solutions are unstable to radiation waves. The region 
of stability of the soliton lies between the straight line and 'parabola', to the right 
of the thick solid curve. The dashed curve gives the variational approximation 
to the stability boundary: h = (1 + 7 4 /4) 1/2 , 7 > a/2. 

where x is the phase of the field ip: ip = \ip\e~ tx . Defining 



whence N(t) — > as t — > 00. We should also mention here that on the other 
side of the wedge, i.e. for h > a/1 + 7 2 , all solutions with ip — > as |x| — > 00 
are unstable against nonlocalised (continuous spectrum) perturbations [TH] . 

In the absence of the damping and driving, i.e. when h = 7 = 0, all lo- 
calised initial conditions in the two-dimensional NLS equation are known to 
either disperse or blow-up in finite time f3J E>1 . Recently it was shown, how- 
ever, that while the soliton ?/>_ remains unstable for all h and 7, the soliton 
ip + stabilises for sufficiently strong damping and driving [10] . (The smallest 
value of 7 for which this soliton can be stable, is 1.006.) The corresponding 




eq.(JU) implies 



d t N < 2(h - i)N, 



stability chart is shown in FigO Our purpose is to explain, in qualitative 
terms, the stabilization mechanism that is at work here. 

3. To this end, we use the variational approach. Equation (JTJ is derivable 
from the stationary action principle with the Lagrangian 

C = e 27 * Re / {iif) t ip* - | V^l 2 - |Vf + \ip\ 4 - H 2 )d*. (5) 
Choosing a bell-shaped trial function jTTj 

V>= y/Ae- ie -( B+i °> 2 , 
with A, B, 9, and o real functions of t, the Lagrangian (jHJ) reduces to 



£ = e 27 * — 
B 



+ ~~~~T~7 + 4 - hcos((f) + 29) cos< 

IB cos^ a> 2 



(6) 



with tan0 — a/B. The corresponding Euler-Lagrange equations are 



A = 8a A - 2-fA + AhA sin(0 + 29) cos - 2M sin[2(0 + 0)] cos 2 0, (7) 
B = 8aB + 2hBsm(4> + 29)cos4>-2hBsm[2^ + 9)]cos 2 (f), (8) 

9 = 1 + 45 - ^A + 2/icos(0 + 2#)cos0- /i cos[2(0 + #)] cos 2 0, (9) 

o- = 4cr 2 — 45 2 + Ai? — 2hB cos(0 + 2$) cos 

+2/i J Bcos[2(0 + #)]cos 2 0. (10) 

The four-dimensional dynamical system defined by (|7|) -(|10 j) has two fixed 
points representing the two stationary solitons: 

A ± = 2(l±^- 7 2 ), B± = ^± 

a 1 • T a 1 • 7 n 
9 + = - arcsm — , 9- = n arcsm — , cr± = 0. 

+ 2 h 2 h 

Consistently with the stability properties of the solitons in the full PDE 
(|T|). the fixed point B_) is unstable for all h and 7 whereas the point 
(A + , B + ) is unstable for small 7 but stabilises for larger dampings. (More 
precisely, this stationary point is stable in the region described by h > 



1 + 7 4 /4, with 7 > y/2 — see Fig^) Therefore, the four-mode approx- 
imation captures the essentials of the infinite-dimensional dynamics in the 
localised- waveform sector. We will now establish two constraints reducing 
the number of independent degrees of freedom to two; these constraints will 
eventually provide the key to the stability mechanism. 
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4. The two-dimensional reduction arises in the overdamped limit, i.e. for 
large 7. In this limit, the dynamics should occur on a slow time-scale; hence 
we introduce the "slow" time T = t/j. We can also expand the solution in 
powers of the small parameter 7 _1 : 

A = A + -A 1 + ..., B = B + -B 1 + 

7 7 

T l n 1 
= - + -01 + ..., <7=-<7 1 + .... 

4 7 7 

Letting /i = 7 + 0/(27) with < c < 1, we make sure that h lies in the 
region of interest: 7 < h < \/l + 7 2 . Substituting in (J7j)- (fTUj) and matching 
coefficients of like powers of 7 -1 , yields a two-dimensional system 

^ = A [c + 8a 1 -4^ + 2 ( ( T 1 / J B ) 2 ], (11) 



where 



8^50 + 4^+4(^/50), (12) 



= ^ + 2B -^A , (13) 
= ^ 5 -25 2 . (14) 



Like their parent system, eqs. (jllj) - (jT2j) have two stationary points in the first 
quadrant of the (A , £? )-plane, 

B o = — Y~i A o = ^ B ti 

with Of = and af = 0. (Since the Aq- and -Bo-axis are invariant 

manifolds, we can restrict our attention to the first quadrant only. In partic- 
ular, fixed points with A < or B < can have no effect on the dynamics 
in the first quadrant.) Like in system (ffj)- (|T0|) with large 7, the fixed point 
(Aq, Bq) is stable (a stable focus) and the other fixed point, (Aq , B^), is 
unstable (a saddle). The sink at the origin is another attractor in the sys- 
tem, competing with the "soliton" (Aq,Bq). The corresponding basins of 
attraction are separated by the stable manifold of the saddle (see Fig|2I). As 
t — > —00, the top part of the separatrix (along with many other trajectories) 
satisfies A , B — > 00 with A oc B^ 2 . 

The most important conclusion of the finite-dimensional analysis is sum- 
marised by eqs. ()13)) - (fT4*J) : two of the four modes are enslaved by the other 
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Figure 2: Phase portrait of the vector field (TTTl)-(TT2l. (In this plot, c = 1/3.) 
The stable manifold of the saddle (Aq,Bq) separates the plane into the basins 
of attraction of the node at the origin and focus at (Aq, B£). 
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two. It is fitting to note here that the division into the "masters" and "slaves" 
is somewhat arbitrary; although in (jli p -([14j ) the amplitude and width ap- 
pear as "master" modes and the two components of the phase as "slaves", 
the reduction can be easily reformulated in such a way that, for example, 9 
and a are the masters and A and B are the slaves. All one needs to do is 
express A and B through 9\ and <j\ from (fT3j) and (HU), and substitute into 
(fTTjl and PI). 

5. In order to explain the stability mechanism, we turn to equation 
P|) governing the density of the soliton's elementary constituents, l^] 2 . [If 
eq.(0) is used to model Faraday resonance in granular media, / IV'I 2 dx has 
the meaning of the total number of particles captured in the oscillon.] The 
first term on the right-hand side of (@J) does not affect the total number 
of the constituents. All it does is rearranges the constituents across the 
oscillon. The second term, on the contrary, does give rise to the creation and 
annihilation of particles. Since this term is proportional to |^| 2 , the creation 
and annihilation occur mainly in the core of the oscillon, where |?/>| 2 is not 
small. In the core we have r ~ and so the creation and annihilation is 
controlled by 9, the uniform component of the phase 

x = 9 + ar 2 . (15) 

The nonuniform part of the phase, or 2 , is small in the core and plays a 
secondary role in this process. Instead, the significance of the quantity a 
is in that it controls the flux of the constituents between the core and the 
periphery of the soliton — see the Vx-term in the r.h.s. of (@J). 

If we perturb the stationary point ip + in the 4-dimensional phase space 
of (fTjl- lfTUj) . the variables 9 and a will zap, within a very short time At ~ -, 
onto the 2-dimensional subspace defined by the constraints fT3|) - fT4|) . After 
this short transient, the evolution of 9 and a will be immediately following 
that of the soliton's amplitude \[A and width 1/^/B. Since the phase x 1S 
coupled to the driver, this provides a negative feedback: perturbations in A 
and B produce only such changes in the two parts of the phase that the new 
values of x{ r ) stimulate the recovery of the stationary values of A and B. 
(The flat phase 9 works to restore the number of constituents while the chirp 
u rearranges them across the soliton.) 

This can be illustrated by considering ip as the envelope of an oscillon 
on the surface of a granular layer, e.g. a layer of tiny brass beads — as in 
the original experiment j2J. Imagine that we increase the amplitude \f~A of 
the oscillon (for example, by dropping several beads on its top): SA (0) > 0. 
Assume, for simplicity, that we do this without changing the oscillon's width: 
8B (0) = 0. From eqs.flESJ,© it follows then that 59 1 (0) < 0. Since 
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the angle 2\ = 2(7r/4 + 7~ 1 # ] f ) is acute for the stationary point (Aq,Bq) 
(remember, 9^ < and erf = 0), the decrease in 9 produces a decrease in 
sin 2%. (Here \ is given by eq. (jl5j) .) As a result, the second term in the 
right-hand side of eq. (jlj) becomes negative which triggers the annihilation of 
elementary constituents. (In the experimental situation this simply means 
that the oscillon starts "leaking" beads to the surrounding medium.) The 
annihilation continues until the original, stationary, value of x (and hence, 
the original value of A = A + ) is recovered. 

Why does this mechanism not work in the case of the unstable fixed 
point, (Aq,Bq)7 The difference is that in that case, 9± > and so 2\ = 
2(7r/4 + 7~ 1 6 l f ) is an obtuse angle. Therefore adding particles at the initial 
moment of time and the resulting decrease of the phase x gi ye r i se to an 
increase of sin 2%. The second term in the right-hand side of (jlj) becomes 
positive and this triggers a further creation of elementary constituents (that 
is, more brass beads will be pulled into the oscillon from the surrounding 
layer.) Hence this time the feedback is positive which makes the stabilisation 
impossible. 

Finally, why is the large damping essential for stability? For small 7 
the coupling of 9 and a to A and B is via differential rather than algebraic 
equations. This time, the dynamics of 9 and a is inertial and so the evolution 
of the phase may not catch up with that of the amplitude and width. The 
feedback loop breaks down and the soliton destabilises. 

In conclusion, the stabilisation mechanism comprises two main ingredi- 
ents: (a) the enslaving of two essential degrees of freedom (e.g. the flat and 
quadratic components of the phase) by another two (amplitude and effec- 
tive width of the soliton); and (b) locking of the phase (and thereby of the 
amplitude and width) to the parametric driver. 
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